Search for: All records

Abstract Large scientific institutions, such as the Space Telescope Science Institute, track the usage of their facilities to understand the needs of the research community. Astrophysicists incorporate facility usage data into their scientific publications, embedding this information in plain text. Traditional automatic search queries prove unreliable for accurate tracking due to the misidentification of facility names in plain text. As automatic search queries fail, researchers are required to manually classify publications for facility usage, which consumes valuable research time. In this work, we introduce a machine learning classification framework for the automatic identification of facility usage of observation sections in astrophysics publications. Our framework identifies sentences containing telescope mission keywords (e.g., Kepler and TESS) in each publication. Subsequently, the identified sentences are transformed using term frequency–inverse document frequency and classified with a support vector machine. The classification framework leverages the context surrounding the identified telescope mission keywords to provide relevant information to the classifier. The framework successfully classifies the usage of MAST-hosted missions with a 92.9% accuracy. Furthermore, our framework demonstrates robustness when compared to other approaches, considering common metrics and computational complexity. The framework’s interpretability makes it adaptable for use across observatories and other scientific facilities worldwide. 

Recent observations have found a growing number of hypervelocity stars with speeds of ≈1500 − 2500 km s⁻¹that could have only been produced through thermonuclear supernovae in white dwarf binaries. Most of the observed hypervelocity runaways in this class display a surprising inflated structure: their current radii are roughly an order of magnitude greater than they would have been as white dwarfs filling their Roche lobe. While many simulations exist studying the dynamical phase leading to supernova detonation in these systems, no detailed calculations of the long-term structure of the runaways have yet been performed. We used an existing AREPOhydrodynamical simulation of a supernova in a white dwarf binary as a starting point for the evolution of these stars with the one-dimensional stellar evolution code MESA. We show that the supernova shock is not energetic enough to inflate the white dwarf over timescales longer than a few thousand years, significantly shorter than the 10^5 − 6year lifetimes inferred for observed hypervelocity runaways. Although they experience a shock from a supernova less than ≈0.02 R_⊙away, our models do not experience significant interior heating, and all contract back to radii of around 0.01 R_⊙within about 10⁴years. Explaining the observed inflated states requires either an additional source of significant heating or some other physics that is not yet accounted for in the subsequent evolution. 

Abstract Stripped-envelope supernovae (SESNe) represent a significant fraction of core-collapse supernovae, arising from massive stars that have shed their hydrogen and, in some cases, helium envelopes. The origins and explosion mechanisms of SESNe remain a topic of active investigation. In this work, we employ radiative-transfer simulations to model the light curves and spectra of a set of explosions of single, solar-metallicity, massive Wolf–Rayet stars with ejecta masses ranging from 4 to 11M_⊙, which were computed from a turbulence-aided and neutrino-driven explosion mechanism. We analyze these synthetic observables to explore the impact of varying ejecta mass and helium content on observable features. We find that the light curve shape of these progenitors with high ejecta masses is consistent with observed SESNe with broad light curves but not the peak luminosities. The commonly used analytic formula based on rising bolometric light curves overestimates the ejecta mass of these high-initial-mass progenitor explosions by a factor of up to 2.6. In contrast, the calibrated method by Haynie et al., which relies on late-time decay tails, reduces uncertainties to an average of 20% within the calibrated ejecta mass range. Spectroscopically, the He i1.083μm line remains prominent even in models with as little as 0.02M_⊙of helium. However, the strength of the optical He ilines is not directly proportional to the helium mass but instead depends on a complex interplay of factors such as the⁵⁶Ni distribution, composition, and radiation field. Thus, producing realistic helium features requires detailed radiative transfer simulations for each new hydrodynamic model. 

Abstract Understanding the evolution of massive binary stars requires accurate estimates of their masses. This understanding is critically important because massive star evolution can potentially lead to gravitational-wave sources such as binary black holes or neutron stars. For Wolf–Rayet (WR) stars with optically thick stellar winds, their masses can only be determined with accurate inclination angle estimates from binary systems which have spectroscopic M sin i measurements. Orbitally phased polarization signals can encode the inclination angle of binary systems, where the WR winds act as scattering regions. We investigated four Wolf–Rayet + O star binary systems, WR 42, WR 79, WR 127, and WR 153, with publicly available phased polarization data to estimate their masses. To avoid the biases present in analytic models of polarization while retaining computational expediency, we used a Monte Carlo radiative-transfer model accurately emulated by a neural network. We used the emulated model to investigate the posterior distribution of the parameters of our four systems. Our mass estimates calculated from the estimated inclination angles put strong constraints on existing mass estimates for three of the systems, and disagree with the existing mass estimates for WR 153. We recommend a concerted effort to obtain polarization observations that can be used to estimate the masses of WR binary systems and increase our understanding of their evolutionary paths. 

Abstract While the Milky Way nuclear star cluster (MW NSC) has been studied extensively, how it formed is uncertain. Studies have shown it contains a solar and supersolar metallicity population that may have formed in situ, along with a subsolar-metallicity population that may have formed via mergers of globular clusters and dwarf galaxies. Stellar abundance measurements are critical to differentiate between formation scenarios. We present new measurements of [M/H] and α -element abundances [ α /Fe] of two subsolar-metallicity stars in the Galactic center. These observations were taken with the adaptive-optics-assisted high-resolution ( R = 24,000) spectrograph NIRSPEC in the K band (1.8–2.6 micron). These are the first α -element abundance measurements of subsolar-metallicity stars in the MW NSC. We measure [M/H] = − 0.59 ± 0.11, [ α /Fe] = 0.05 ± 0.15 and [M/H] = − 0.81 ± 0.12, [ α /Fe] = 0.15 ± 0.16 for the two stars at the Galactic center; the uncertainties are dominated by systematic uncertainties in the spectral templates. The stars have an [ α /Fe] in between the [ α /Fe] of globular clusters and dwarf galaxies at similar [M/H] values. Their abundances are very different than the bulk of the stars in the nuclear star cluster. These results indicate that the subsolar-metallicity population in the MW NSC likely originated from infalling dwarf galaxies or globular clusters and are unlikely to have formed in situ. 

Abstract Type Ia supernovae (SNe Ia) are securely understood to come from the thermonuclear explosion of a white dwarf as a result of binary interaction, but the nature of that binary interaction and the secondary object is uncertain. Recently, a double white dwarf model known as the dynamically driven double-degenerate double-detonation (D6) model has become a promising explanation for these events. One realization of this scenario predicts that the companion may survive the explosion and reside within the remnant as a fast moving (V_peculiar> 1000 km s⁻¹), overluminous (L> 0.1L_⊙) white dwarf. Recently, three objects that appear to have these unusual properties have been discovered in the Gaia survey. We obtained photometric observations of the SN Ia remnant SN 1006 with the Dark Energy Camera over four years to attempt to discover a similar star. We present a deep, high-precision astrometric proper-motion survey of the interior stellar population of the remnant. We rule out the existence of a high-proper-motion object consistent with our tested realization of the D6 scenario (V_transverse> 600 km s⁻¹withm_r< 21 corresponding to an intrinsic luminosity ofL> 0.0176L_⊙). We conclude that such a star does not exist within the remnant or is hidden from detection by either strong localized dust or the unlikely possibility of ejection from the binary system almost parallel to the line of sight. 

Abstract With the advent of high-cadence, all-sky automated surveys, supernovae (SNe) are now discovered closer than ever to their dates of explosion. However, young premaximum light follow-up spectra of Type Ic SNe (SNe Ic), probably arising from the most-stripped massive stars, remain rare despite their importance. In this Letter, we present a set of 49 optical spectra observed with the Las Cumbres Observatory through the Global Supernova Project for 6 SNe Ic, including a total of 17 premaximum spectra, of which 8 are observed more than a week before V -band maximum light. This data set increases the total number of publicly available premaximum-light SN Ic spectra by 25%, and we provide publicly available SNID templates that will significantly aid in the fast identification of young SNe Ic in the future. We present a detailed analysis of these spectra, including Fe ii 5169 velocity measurements, O i 7774 line strengths, and continuum shapes. We compare our results to published samples of stripped SNe in the literature and find one SN in our sample that stands out. SN 2019ewu has a unique combination of features for an SN Ic: an extremely blue continuum, high absorption velocities, a P Cygni–shaped feature almost 2 weeks before maximum light that TARDIS radiative transfer modeling attributes to C ii rather than H α , and weak or nonexistent O i 7774 absorption feature until maximum light. 

Abstract Type Iax supernovae (SNe Iax) are the largest known class of peculiar white dwarf SNe, distinct from normal Type Ia supernovae (SNe Ia). The unique properties of SNe Iax, especially their strong photospheric lines out to extremely late times, allow us to model their optical spectra and derive the physical parameters of the long-lasting photosphere. We present an extensive spectral timeseries, including 21 new spectra, of SN Iax 2014dt from +11 to +562 days after maximum light. We are able to reproduce the entire timeseries with a self-consistent, nearly unaltered deflagration explosion model from Fink et al. using TARDIS , an open source radiative-transfer code. We find that the photospheric velocity of SN 2014dt slows its evolution between +64 and +148 days, which closely overlaps the phase when we see SN 2014dt diverge from the normal spectral evolution of SNe Ia (+90 to +150 days). The photospheric velocity at these epochs, ∼400–1000 km s −1 , may demarcate a boundary within the ejecta below which the physics of SNe Iax and normal SNe Ia differ. Our results suggest that SN 2014dt is consistent with a weak deflagration explosion model that leaves behind a bound remnant and drives an optically thick, quasi-steady-state wind creating the photospheric lines at late times. The data also suggest that this wind may weaken at epochs past +450 days, perhaps indicating a radioactive power source that has decayed away. 

Abstract We study the production of very light elements ( Z < 20) in the dynamical and spiral-wave wind ejecta of binary neutron star mergers by combining detailed nucleosynthesis calculations with the outcome of numerical relativity merger simulations. All our models are targeted to GW170817 and include neutrino radiation. We explore different finite-temperature, composition-dependent nuclear equations of state, and binary mass ratios, and find that hydrogen and helium are the most abundant light elements. For both elements, the decay of free neutrons is the driving nuclear reaction. In particular, ∼0.5–2 × 10 −6 M ⊙ of hydrogen are produced in the fast expanding tail of the dynamical ejecta, while ∼1.5–11 × 10 −6 M ⊙ of helium are synthesized in the bulk of the dynamical ejecta, usually in association with heavy r -process elements. By computing synthetic spectra, we find that the possibility of detecting hydrogen and helium features in kilonova spectra is very unlikely for fiducial masses and luminosities, even when including nonlocal thermodynamic equilibrium effects. The latter could be crucial to observe helium lines a few days after merger for faint kilonovae or for luminous kilonovae ejecting large masses of helium. Finally, we compute the amount of strontium synthesized in the dynamical and spiral-wave wind ejecta, and find that it is consistent with (or even larger than, in the case of a long-lived remnant) the one required to explain early spectral features in the kilonova of GW170817. 

We present the first results of a comprehensive supernova (SN) radiative-transfer (RT) code-comparison initiative (StaNdaRT), where the emission from the same set of standardised test models is simulated by currently used RT codes. We ran a total of ten codes on a set of four benchmark ejecta models of Type Ia SNe. We consider two sub-Chandrasekhar-mass (M_tot= 1.0M_⊙) toy models with analytic density and composition profiles and two Chandrasekhar-mass delayed-detonation models that are outcomes of hydrodynamical simulations. We adopt spherical symmetry for all four models. The results of the different codes, including the light curves, spectra, and the evolution of several physical properties as a function of radius and time are provided in electronic form in a standard format via a public repository. We also include the detailed test model profiles and several Python scripts for accessing and presenting the input and output files. We also provide the code used to generate the toy models studied here. In this paper, we describe the test models, radiative-transfer codes, and output formats in detail, and provide access to the repository. We present example results of several key diagnostic features.